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5 CROSS REFERENCE TO RELATED APPLICATIONS 
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and U.S. Provisional Application 60/462,314, filed 14 
April 2003. The contents of these applications are 
10 incorporated herein by reference. 

STATEMENT ON FEDERALLY SPONSORED RESEARCH 

N/A 

15 FIELD OF THE INVENTION 

The present invention relates to a mass 
spectrometer and a method of mass spectrometry. 

BACKGROUND OF THE INVENTION 
2 0 Known mass spectrometers measure and report a value 

of mass to charge ratio for analyte ions in a sample. 
The determined mass to charge ratio is taken to be 
representative of the sum of the masses of the 
individual atoms present in an analyte molecule plus or 

2 5 minus the rest mass of one or more electrons (depending 

upon the charge state and polarity of the molecule when 
measured) . The unit of mass to charge is based upon one 
mass unit being 1/12 of the mass of the most abundant 
isotope of carbon ( 12 C) and one unit of charge being 

3 0 equal to the charge of an electron. 

Conventional mass spectrometers process mass 
spectral data and calculate the mass to charge ratio of 
detected ions. However, the accuracy of a mass to 
charge ratio measurement will be dependent upon the 
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particular type of mass analyser used, the quantity of 
analyte sample and the conditions under which a 
particular sample is mass analysed. 

Although some conventional mass spectrometers are 
5 able to measure the mass to charge ratio of an ion to a 
relatively high precision and accuracy, conventional 
mass spectrometers are not s able to determine the 
accuracy of each individual mass to charge ratio 
measurement relating to an unknown sample. 

10 However, without an estimate of the accuracy of a 

mass to charge ratio measurement it is uncertain what 
range of mass to charge ratio should be considered when 
using, for example, a database of atomic elements and 
their isotopes to determine possible elemental 

15 compositions which will have mass to charge ratios 
substantially similar to that of the analyte ion. 
Similarly, it will also not be known what range of mass 
to charge ratio should be considered when searching a 
database of known molecules to find a match for an 

2 0 analyte ion. It will be appreciated that if a 

relatively wide mass to charge ratio window needs to be 
considered then the number of possible candidate ions 
dramatically increases . 

It is therefore desired to provide an improved 

2 5 method of mass spectrometry and an improved mass 

spectrometer . 

SUMMARY OF THE INVENTION 

3 0 According to an aspect of the present invention 

there is provided a mass spectrometer comprising: 
a mass analyser; and 

a processing system, the processing system being 
arranged to obtain mass spectral data, to determine a 
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mass to charge ratio of n different species of ion 
observed in the mass spectral data and to calculate a 
respective error band for the determined mass to charge 
ratio of each of the n different species of ion, wherein 
5 the probability or confidence that the real, true, 

actual or accepted mass to charge ratio of a species of 
ion falls within its respective calculated error band is 
greater than or equal to x% . 

According to different embodiments of the present 

10 invention n may be 1, 2-5, 5-10, 10-15, 15-20, 20-25, 

25-30, 30-35, 35-40, 40-45, 45-50, > 50, > 60, ^ 70, > 
80, > 90 or > 100. 

According to less preferred embodiments x may be 
<1, 1-5, 5-10, 10-15, 15-20, 20-25, 25-30, 30-35, 35-40, 

15 40-45, 45-50. According to more preferred embodiments x 
may be 50-55, 55-60, 60-65, 65-70, 70-75, 75-80, 80-85 
or 85-90. According to particularly preferred 
embodiments x may be 90-91, 91-92, 92-93, 93-94, 94-95, 
95-96, 96-97, 97-98, 98-99, 99-99.5, 99.5-99.95, 99.95- 

20 99.99 or 99.99-100. 

The mass spectrometer may further comprise means 
for, reporting the mass to charge ratio of at least' some 
of the n species of ion together with the respective 
calculated error band for each mass to charge ratio. 

25 The mass spectral data preferably comprises data 

relating to the time that different species of ion were 
detected. The data relating to the time that different 
species of ion were detected is preferably converted 
into data relating to the mass to charge ratio of the 

30 ions. The data relating to the time that different 
species of ion were detected is preferably converted 
into data relating to the mass to charge ratio of the 
ions by using a calibration function. 
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As part of the process of calculating a respective 
error band for the determined mass to charge ratio of 
each of the n different species of ion the processing 
system may estimate, in use, the systematic error in 
5 determining the mass to charge ratio of the ions. 

The processing system preferably estimates, in use, 
the accuracy of a mathematical model used in mass 
calibration. According to an embodiment the error due 
to mass calibration may be estimated. 

10 The error due to using one or more internal 

references or calibrants may be estimated. Preferably, 
as part of the process of calculating a respective error 
band for the determined mass to charge ratio of each of 
the n different species of ion the processing system 

15 monitors for changes in the time of detection of one or 
more internal references or calibrants. 

According to an embodiment the stability or 
instability of the mass spectrometer subsequent to mass 
calibration may be estimated. Preferably, as part of 

20 the process of calculating a respective error band for 
the determined mass to charge ratio of each of the n 
different species of ion the processing system 
estimates, in use, the error due to drift of the mass 
spectrometer subsequent to mass calibration. The effect 

25 of changes in the operating conditions of the mass 
spectrometer subsequent to mass calibration may be 
estimated. 

According to another embodiment the effect of user 
intervention subsequent to mass calibration is 
30 estimated. The effect of external influences on the 
mass spectrometer subsequent to mass calibration may 
also be estimated. According to a yet further 
embodiment the effect of the stability or instability of 
a power supply of the mass spectrometer subsequent to 



5 

Docket No. AE-MM-11 
(MST-1237 US2) 

mass calibration is estimated. Preferably, as part of 
the process of calculating a respective error band for 
the determined mass to charge ratio of each of the n 
different species of ion the processing system 
5 estimates, in use, the effect of changes in temperature 
on the mass spectrometer subsequent to mass calibration. 

The effect of changes in interference subsequent to 
mass calibration may be estimated. The interference may 
be due to one or more internal reference or calibrant 

10 ions interfering with one or more analyte ions. The 
interference may be due to one or more analyte ions 
interfering with one or more internal reference or 
calibrant ions. The interference may be due to one or 
more background or chemical noise ions interfering with 

15 one or more internal reference or calibrant ions. 

Alternatively/additionally, the interference may be due 
to one or more background or chemical noise ions 
interfering with one or more analyte ions. 

As part of the process of calculating a respective 

2 0 error band for the determined mass to charge ratio of 

each of the n different species of ion the processing 
system may estimate, in use, the time elapsed since the 
mass spectrometer was last calibrated. The processing 
system may also estimate, in use, the effect of thermal 
25 expansion of one or more components of the mass 
spectrometer subsequent to mass calibration. 

According to another embodiment the processing 
system may estimate, in use, the error due to space 
charge repulsion effects within the mass spectrometer. 

3 0 The error due to space charge repulsion effects may be 

determined at least in part by consideration of the 
observed intensity of one or more species of ions. 

As part of the process of calculating a respective 
error band for the determined mass to charge ratio of 
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each of the n different species of ion the processing 
system may estimate, in use, the error recorded during 
mass calibration. 

According to a particularly preferred embodiment 
the processing system may estimate, in use, the 
statistical error or random error in determining the 
mass to charge ratio of the ions. In particular, the 
processing system may estimate, in use, the uncertainty 
or standard deviation of the time of detection of one or 
more species of ion. According to a less preferred 
embodiment, particularly using a Fourier Transform Ion 
Cyclotron Resonance mass analyser the processing system 
may estimate, in use, the uncertainty or standard 
deviation of the frequency of detection of one or more 
species of ion. 

Preferably, as part of the process of calculating a 
respective error band for the determined mass to charge 
ratio of each of the n different species of ion the 
processing system estimates, in use, the error due to 
ion detection statistics. The processing system may 
estimate, in use, the error due to insufficient 
sampling . 

According to another embodiment, as part of the 
process of calculating a respective error band for the 
determined mass to charge ratio of each of the n 
different species of ion the processing system 
estimates, in use, the error due to computational 
errors. Preferably, as part of the process of 
calculating a respective error band for the determined 
mass to charge ratio of each of the n different species 
of ion the .processing system estimates, in use, the 
error due to rounding off errors. 

The processing system may estimate, in use, the 
error in a report exception situation wherein one or 
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more species of ion within the mass spectral data are 
sufficiently corrupted or suffering from interference 
that the calculated respective error band of the ions is 
not reported. 

5 As part of the process of calculating a respective 

error band for the determined mass to charge ratio of 
each of the n different species of ion the processing 
system may estimate, in use, the error due to mass 
interference wherein background or interference ions 

10 having substantially similar mass to charge ratios to 

analyte ions affect the ability of resolving the analyte 
ions from the background or interference ions . The 
error due to mass interference may be estimated by 
comparing the shape of an observed mass peak from the 

15 mass spectral data with the shape of a theoretical mass 
peak. The error due to mass interference may, for 
example, be estimated by comparing the deviation of the 
shape of an observed mass peak with either a Gaussian 
distribution, a cosine squared distribution, another 

20 mathematical function, or a profile derived from the 

shape of one or more mass peaks of one or more internal 
reference or calibrant ions. 

According to another embodiment as part of the 
process of calculating a respective error band for the 

2 5 determined mass to charge ratio of each of the n 
different species of ion the processing system 
estimates, in use, the error due to saturation of a 
detector. According to another embodiment the 
processing system may estimate, in use, the error due to 

30 signal amplification saturation. According to another 
embodiment the processing system may estimate, in use, 
the error due to deadtime effects or missed counts of an 
ion counting detector. 
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As part of the process of calculating a respective 
error band for the determined mass to charge ratio of 
each of the n different species of ion the processing 
system may estimate, in use, the error due to the 
5 response of an ion detector. 

According to another embodiment the processing 
system estimates, in use, the error due to the response 
of an electron or photon multiplier. 

The processing system preferably combines multiple 

10 estimates of different errors. The multiple estimates 
of different errors are preferably combined by adding 
the estimates in quadrature. 

According to a preferred embodiment the mass 
spectrometer further comprises means arranged to search 

15 in use a database to look for one or more ions known to 
have mass to charge ratios within a range which is 
related to the determined mass to charge ratio of one or 
more ions and the respective calculated error band for 
the determined mass to charge ratio. The database may 

20 comprise details of biopolymers, proteins, peptides, 
polypeptides , oligionucleotides , oligionucleosides , 
amino acids, carbohydrates, sugars, lipids, fatty acids, 
vitamins, hormones , portions or fragments of DNA, 
portions or fragments of cDNA, portions or fragments of 

25 RNA, portions or fragments of mRNA, portions or 

fragments of tRNA, polyclonal antibodies, monoclonal 
antibodies, ribonucleases , enzymes, metabolites, 
polysaccharides, phosphorolated peptides, phosphorolated 
proteins, glycopeptides , glycoproteins or steroids. 

30 Additionally/alternatively, the database may comprise 
details of the Electron Impact mass spectra of 
compounds. For example, the Wiley and/or NBS libraries 
may be searched. 
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According to an embodiment the mass spectrometer 
further comprises means for calculating elemental 
compositions having mass to charge ratios within a 
calculated error band of the determined mass to charge 
5 ratio of one or more species of ion. The elemental 

compositions are preferably calculated using the known 
masses or mass to charge ratios of atomic elements or 
groups of elements and/or their isotopes. 

According to an embodiment the processing system is 

10 arranged, in use, to define one or more compounds of 
interest and to interrogate at least some of the 
determined mass to charge ratios and respective 
calculated error bands of at least some of the different 
species of ion observed in the mass spectral data to see 

15 whether the one or more compounds of interest fall 

within the error band of at least some of the different 
species of ion. 

According to an embodiment the processing system is 
arranged, in use, to define one or more mass to charge 

20 ratios of interest and to interrogate at least some of 
the determined mass to charge ratios and respective 
calculated error bands of at least some of the different 
species of ion observed in the mass spectral data to see 
whether the one or more mass to charge ratios of 

25 interest fall within the error band of at least some of 
the different species of ion. 

According to an embodiment the processing system is 
arranged, in use, to define one or more compounds of 
interest and to interrogate at least some of the 

30 determined mass to charge ratios and respective 

calculated error bands of different species of ion 
observed from a plurality of mass spectra acquired over 
a period of time to see whether the one or more 
compounds of interest fall within the error band of at 
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least some of the different species of ion observed from 
the plurality of mass spectra. 

According to an embodiment the processing system is 
arranged, in use, to define one or more mass to charge 
5 ratios of interest and to interrogate at least some of 
the determined mass to charge ratios and respective 
calculated error bands of different species of ion 
observed from a plurality of mass spectra acquired over 
a period of time to see whether the one or more mass to 

10 charge ratios of interest fall within the error band of 
at least some of the different species of ion observed 
from the plurality of mass spectra. 

The mass spectrometer preferably further comprises 
an Electrospray ("ESI") ion source, an Atmospheric 

15 Pressure Chemical Ionisation ("APCI") ion source, an 
Atmospheric Pressure Photo Ionisation ("APPI") ion 
source, a Laser Desorption Ionisation ( "LDI " ) ion 
source, an Inductively Coupled Plasma ("ICP") ion - 
source, an Electron Impact ("EI) ion source, a Chemical 

20 Ionisation ("CI") ion source, a Field Ionisation ("FI") 
ion source, a Fast Atom Bombardment ("FAB") ion source, 
a Liquid Secondary Ion Mass Spectrometry ("LSIMS" ) ion 
source, an Atmospheric Pressure Ionisation ("API") ion 
source, a Field Desorption ("FD") ion source, a Matrix 

25 Assisted Laser Desorption Ionisation ("MALDI") ion 

source or a Desorption/Ionisation on Silicon ("DIOS") 
ion source. The ion source may be a continuous ion 
source or a pulsed ion source. 

The mass spectrometer preferably comprises a Time 

30 of Flight mass analyser, a quadrupole mass analyser, a 
Penning or Fourier Transform Ion Cyclotron Resonance 
("FTICR") mass analyser, a 2D or linear quadrupole ion 
trap or a Paul or 3D quadrupole ion trap. 
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According to another aspect of the present 
invention there is provided a method of mass 
spectrometry comprising: 

obtaining mass spectral data using a mass 
5 spectrometer; 

determining a mass to charge ratio of n different 
species of ion observed in the mass spectral data; and 

calculating a respective error band for the 
determined mass to charge ratio of each of the n 
10 different species of ion, wherein the probability or 

confidence that the real, true, actual or accepted mass 
to charge ratio of a species of ion falls within its 
respective calculated error band is greater than or 
equal to x% . 

15 The preferred embodiment is particularly 

advantageous in that an understanding of the error in a 
mass to charge ratio measurement and the overall 
confidence in the mass to charge ratio assignment for an 
ion is particularly useful if the mass assignment is to 

2 0 be used, for example, to infer information about the 

elemental composition of the ion. 

The preferred embodiment relates to automatically 
estimating and optionally reporting the expected or 
estimated accuracy of individual mass to charge ratio 
25 measurements. This is not something which conventional 
mass spectrometers are capable of doing. 

The preferred embodiment automatically provides an 
estimate of the accuracy to be expected for individual 
mass to charge ratio measurements of unknown samples. 

3 0 This provides accurate limits of mass to charge ratio 

within which candidate elemental compositions may lie 
in, for example, a database of analyte ions. Estimating 
the accuracy of an assignment of mass to charge ratio to 
observed analyte ions also provides a measure of 
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confidence in such elemental composition assignments. 
It provides both a means of restricting the range of 
possible elemental compositions which may need to be 
considered and also guards against over restriction 
5 thereby reducing the possibility of erroneous assignment 
of elemental composition. The preferred embodiment is 
also particularly advantageous when searching against 
mass to charge ratios in a library of known ions having 
known mass to charge ratios. 

10 The precision of a mass to charge ratio measurement 

is dependent upon the scatter of measured values and the 
number of measurements. For a normal distribution of 
measurements the precision may be estimated from the 
standard deviation a of measurements about the mean and 

15 the number of measurements n. The accuracy of a 

particular mass to charge ratio measurement is the 
difference between the mean of the mass to charge ratio 
measurements and the true, actual, real or accepted mass 
to charge ratio of the analyte ion. In general, the 

2 0 more accurate the measurement of mass to charge ratio 

the smaller the band of uncertainty and hence the more 
restricted the range of possible elemental compositions 
which have mass to charge ratios which fall within this 
band of uncertainty. 
25 Conventional mass spectrometers and in particular 

Fourier Transform Ion Cyclotron Resonance mass 
spectrometers, double focusing magnetic sector mass 
spectrometers, orthogonal and axial Time of Flight mass 
spectrometers and RF quadrupole mass filter mass 

3 0 spectrometers are capable of measuring the mass to 

charge ratio of an ion to within a few parts per million 
(ppm) of the true mass to charge ratio. Although such 
mass spectrometers are able to measure mass to charge 
ratio with a high degree of precision and accuracy, no 
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estimate of the accuracy of individual mass to charge 
ratio measurements is possible. 

Without an estimate of the accuracy of a mass to 
charge ratio measurement the value or worth of such a 
5 mass to charge ratio measurement is unknown. An 

understanding of the possible error in a mass to charge 
ratio measurement and an estimate of the confidence or 
otherwise in the mass to charge ratio assignment for an 
individual ion is particularly beneficial if the mass to 

10 charge ratio assignment is to be used to infer 

information about the elemental composition of the ion. 

The preferred embodiment allows the mass to charge 
ratio window used to calculate proposed elemental 
compositions to be automatically restricted. 

15 The preferred estimation of mass to charge ratio 

accuracy may also be used to appropriately set a search 
window when trying to match the measured mass to charge 
ratio or mass of an analyte ion to details of ions or 
molecules held in a. library of ions or molecules having 

20 known mass to charge ratios or masses. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Various embodiments of the present invention will 
2 5 now be described, by way of example only, and with 
reference to the accompanying drawings in which: 

Fig. 1 shows a theoretical mass spectral peak; 

Fig. 2 shows the theoretical mass spectral peak in 
which the mass to charge ratio axis has been divided 
30 into discrete samples; 

Fig. 3 shows an uncalibrated mass spectrum of a 
reference compound PFTBA; 

Fig. 4 shows how residual errors after mass 
calibration may be quite high if the mass spectrometer 
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is only calibrated using a relatively small amount of 
data; 

Fig. 5 shows how the residual errors after mass 
calibration will be smaller but not insignificant if the 
5 mass spectrometer is calibrated using a relatively large 
amount of data; 

Fig. 6 shows a portion of a mass spectrum wherein a 
resolution of 10,000 is sufficient to resolve two 
different ions; 

10 Fig. 7 shows a portion of a mass spectrum wherein 

two different ions can not be resolved; 

Fig. 8 shows the residual errors and calculated 95% 

confidence interval after calibration according to the 

preferred embodiment; and 
15 Fig. 9 shows a mass spectrum according to the 

preferred embodiment wherein the mass to charge ratio of 

individual ions together with the respective error in 

the determination of the mass to charge ratio of the ion 

is reported. 

20 

DETAILED DESCRIPTION OF THE INVENTION 

A preferred embodiment of the present invention 
will now be described. The preferred embodiment 

2 5 involves arriving at a representative value for the 

precision and accuracy of an individual mass to charge 
ratio measurement. In order to do this multiple sources 
of error are preferably estimated and at least some of 
these errors are then combined and reported along with 

30 each individual mass to charge ratio measurement. 

Some of the sources of error involved with an 
individual mass to charge ratio measurement which are 
preferably estimated include statistical or random 
errors such as errors due to ion detection statistics, 
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errors due to insufficient sampling frequency, 
computational errors (e.g. rounding off errors) and 
errors recorded during mass calibration. 

In some cases data may be sufficiently corrupted 
5 that no estimate of mass measurement error may be 
reliably calculated or reported. These situations 
include the effects of mass interferences, detector 
saturation and errors due to instability or drift of the 
mass spectrometer subsequent to calibration or in 

10 between calibrations. 

Random error in the determination of the mass to 
charge ratio of an ion is often dominated by ion 
statistics. A mass spectral peak in a mass spectrum is 
the result of detecting a number of ions generated from 

15 an ion source. An example of a theoretical mass peak is 
shown in Fig. 1. Fig. 1 shows a theoretical mass peak 
which may be observed for ions having a mass to charge 
ratio of 300. The detected ions which are all the same 
species of ion will have a distribution of mass to 

2 0 charge ratios centred around a mean value and will fall 
within a defined mass envelope. The narrower the mass 
or peak envelope the more able a mass spectrometer is to 
distinguish between ions with very small differences in 
mass to charge ratio. The mass resolution R may be 

25 defined as: 

SM 

where M is the mass to charge ratio of the mean of a 
30 mass or peak envelope and 5M is the width of the mass or 
peak envelope in mass units defined in a specific way 
e.g. at the Full Width Half Maximum (FWHM) . The mean of 
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this distribution is the estimation of the mass to 
charge ratio for the observed mass peak. 

If the gain characteristics of the detection and 
amplification system are known, then the area of the 
5 mass spectral peak may be expressed in terms of the 
absolute number of ions within the mass or peak 
envelope . 

In Fig. 2 the mass to charge ratio scale has been 
divided into a plurality of discrete samples. Each 

10 sample is described by an index number i and each index 
number i is associated with a value of mass to charge 
ratio Mi. For simplicity, it has been assumed that the 
y axis represents the absolute number of ions Ii 
measured at each particular value of mass to charge 

15 ratio Mi. It is also assumed that each ion detected 

gives rise to a signal of the same amplitude. The mean 
of this distribution M mean is given by: 



M = — 

mean j 

20 

where 



2,1,-* 

i 

2 5 and N is the total number of ions detected. The mean 

represents the measured mass to charge ratio value of 
the peak. 

The mass spectral peak shown in Fig. 2 may be 
considered as a probability distribution with a parent 

3 0 population variance V of ion arrivals which can be 

estimated from: 
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y i;(M:-M m v 

y / j j i v i mean J 

N-l 



The standard deviation a M around the mean for this 
5 parent population is given by: 

From theoretical statistics, if the mass to charge 
10 ratio value is measured based upon a single recorded ion 
belonging to this distribution then the value Mi recorded 
for this ion will be the most likely estimate of the 
mean of the distribution. The estimated uncertainty in 
this single measurement is given by a M . 
15 If the measurement is repeated for N single ions an 

average mean value can be produced. This average mean 
value will have an uncertainty or standard deviation a c 
of: 



20 




where a c represents the standard deviation on the 
measurement of the center of a mass spectral peak with a 
parent population standard deviation of <r M and 

2 5 containing N ions. 

The expressions above can be related to any peak 
shape. However, the shape of a typical mass spectral 
peak may be approximated by a Gaussian distribution. 
Assuming a Gaussian distribution and 8M defined at Full 

3 0 Width Half Maximum (FWHM) then: 
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SM 

cr M = 

M 2.35 

This allows the uncertainty in the centroid 
determination of the peak to be expressed in terms of 5M 
and the number of ions recorded N. The expression below 
yields the uncertainty in the measurement in mass units: 

SM 

From integration of the Gaussian probability 
distribution function, approximately 95% of measurements 
for a given SM and N will fall within +/- 2a c . This 
spread is referred to as the 95% confidence limit. 

The equations derived above describe the 
statistical limitations of peak centroid determination 
in a situation where every ion detected gives an 
identical response . 

Typically mass spectrometers employ multiplier 
systems to amplify the signal from the individual ions 
to ensure detection. These multipliers include electron 
multipliers, photo multipliers and microchannel plates 
in a variety of configurations and are often used in 
conjunction with electronic amplification devices. Due 
to the statistical characteristics of some of these 
devices, ions having the same mass to charge ratio will 
have a response governed by a separate probability 
distribution function. When employing an analogue, 
signal recording system (e.g. an Analogue to Digital 
Recorder ADC) this distribution leads to an increase in 
the uncertainty of the determination of the centre of a 
mass spectral peak. This increase may be calculated if 
the distribution of ion response is well characterized 
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and may then be incorporated into the estimation of 
error. When using an ion counting system (e.g. a Time 
to Digital Recorder TDC) this effect does not need to be 
considered. 

5 In some cases there may be insufficient samples 

over a mass spectral peak to accurately determine its 
center. This error should ideally be minimized. 
However, if it is significant then it should preferably 
be incorporated into any final quantification of overall 
10 error. 

When measuring mass to less than 1/1000 of a mass 
unit (commonly referred to as milli dalton (mDa) ) it is 
important that calculations are performed and presented 
with a suitable number of decimal places to avoid 
15 unnecessary rounding errors. For example, if an ion 

having a mass to charge ratio of 50 is mass measured and 
the result rounded to four decimal places then the 
possible error due to this rounding will be 50.0000 ± 
0.00005. Expressed in parts per million (ppm) this is 

2 0 50 ± 1 ppm. 

Calibration involves introducing a reference 
compound, or mixture of reference compounds which give 
rise to mass spectral peaks which have precisely known 
elemental compositions and hence which have precisely 
25 known mass to charge ratios. The experimentally 

measured mass to charge ratio value for each of the 
calibration peaks are then compared to the known 
reference mass to charge ratio values and a suitable 
calibration function or functions can then be applied to 

3 0 subsequent experimental data. The calibration function 

ensures that an unknown sample is correctly calibrated 
as far as is possible. The purpose of this calibration 
function is to adjust the measured mass to charge ratio 
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value of reference ions to be as close as possible to 
the theoretical mass of the reference compound. 

As long as statistical errors are low this 
procedure will compensate for any systematic errors 
5 which are present. Systematic errors in mass 

measurement are biases in the measured value which 
remain the same no matter how many times the measurement 
is repeated. Any residual errors remaining after 
calibration may be treated as random errors. 

10 Ideally, during calibration, sufficient signal 

should be acquired such that statistical error due to an 
insufficient number of ions is negligible. This may 
not, however, always be possible. 

From the statistical expressions developed above it 

15 is possible to assign a value of centroid standard 

deviation to each peak in the calibration. This may be 
used to apply appropriate weighting to each mass to 
charge ratio value when fitting a calibration function 
or functions to the calibration data. This information 

20 may be used to predict how errors due to statistics in 
the initial calibration will affect the uncertainty in 
the final analyte mass measurement. This is illustrated 
further with reference to Figs. 3-5. 

Fig. 3 shows an un-calibrated fragment ion mass 

25 spectrum over the mass to charge ratio range 300-620 
resulting from a reference compound 
Perf luorotributylamine (PFTBA) acquired using an 
Electron Impact ("EI") ion source operated in positive 
ion mode. The ions were analysed using an orthogonal 

30 acceleration Time of Flight mass analyser. PFTBA is a 
common reference material for Electron Impact mass 
spectrometry studies. As a result of the ionisation 
process the molecular ion having a nominal mass to 
charge ratio of 671 was fragmented into a number of 
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different fragment ions having different mass to charge 
ratios. The molecular ion is not therefore observed in 
the resulting fragmentation mass spectrum a portion of 
which is shown in Fig. 3. 
5 With orthogonal acceleration Time of Flight mass 

analysers the mass to charge ratio value is proportional 
to the square of the time of flight of ions in a drift 
or flight region of the mass analyser. By comparing the 
square root of the flight time of a reference ion to the 

10 known mass to charge ratio of the reference ion, a 

relationship between time of flight and mass to charge 
ratio can be established. This relationship may then be 
used to generate a calibration expression based on a 
third order polynomial expression using a weighted least 

15 squares method. 

Fig. 4 shows the residual errors after generating 
such a calibration expression based upon a relatively 
small amount of data. In Fig. 4 the x axis indicates 
the calculated mass to charge ratios of different 

2 0 fragment ions derived from PFTBA based upon their known 
empirical formula. The y axis represents the difference 
between the calibrated mass to charge ratio of the 
fragment ions and the known reference mass to charge 
ratio values after calibration had been applied. The 

25 axis is displayed in parts per million error (ppm) . The 

vertical lines represent +/- 2a c where a c is the 
statistical uncertainty in the peak center determination 
in ppm. These error bars represent the calculated 
scatter which could be expected from the reported mass 
30 to charge ratio value from repetitive measurements of 
each calibration point within a 95% confidence limit. 
The error bars are inversely proportional to the square 
root of the intensity. 
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Fig. 5 shows a similar plot but obtained using a 
much larger data set. The error bars are very small but 
they are nonetheless still present and the total error 
is now less than +/- Ippm over the entire range. The 
5 deviation of the calibrated results from the calculated 
result reflects a remaining element of systematic error. 
This error may be further reduced by employing a higher 
order polynomial expression or an alternative curve 
fitting method. 

10 The results shown in Fig. 4 and 5 were taken from 

the same experimental run. PFTBA was continuously 
infused into the ion source and mass spectra were 
acquired at a rate of one mass spectrum per second. 
Data from multiple mass spectra were averaged prior to 

15 calibration. Fig. 4 shows the error for an average of 
20s of data collection whereas Fig. 5 shows the error 
for an average of 3 0 minutes of data collection. The 
results shown in Fig. 5 therefore represent a data set 
which is approximately x90 that of the data set used in 

2 0 Fig. 4. The error bars, which are based on ion 

statistics alone are approximately x9 . 5 (square root of 
90) times smaller in Fig. 5 than in Fig. 4. 

Because of the variation in signal intensity it can 
be seen from Fig. 4 that the error due to statistics 

25 within the calibration varies for different mass 

regions. The errors in this case are dominated by the 
calculated statistical error on each data point. 

In the specific example of calibrating a Time of 
Flight mass spectrometer coefficients have been 

30 calculated of a third order polynomial of the form: 
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where M re f(i) is the calculated mass to charge ratio 
value of reference peak i based on the known empirical 
formula and M es t(i) is the mass to charge ratio value of 
reference peak i recorded before calibration. 

The method involves minimizing the sum of the 
squared difference between the estimated mass to charge 
ratio values M es t(i) and the calculated mass to charge 
ratio value of the reference material M re f(i> weighted by 
the calculated statistical variance a c( i) 2 of each 
estimated measurement where i is an index denoting 
individual points in the calibration and N is the total 
number of points in the calibration. This sum is 
denoted by the symbol % 2 : 

M refV) 

°"c(0 

where : 

/w« (0 ) = a + b - V M "<«> + c • W M « W ) 2 + d ■ W M «<(o> 3 

20 

The coefficients of the polynomial function are 
adjusted to give the lowest value of % 2 . The % 2 value 
gives an overall estimation of the goodness of fit of 
the calibration expression and can be used to provide 
25 crude acceptance or rejection criteria for a given 
calibration . 

Other goodness of fit tests may be performed on 
this type of data depending on the number of calibration 
points and the calibration approach followed. In 
3 0 principle other methods may be used for rejection or 
acceptance of calibration. 



5 
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Using the method of least squares fitting above, 
the variance in the calculated coefficients in the 
polynomial expression can be calculated along with the 
covariance. This information can be used to calculate 
5 the variance and hence standard deviation of a 

measurement of any mass to charge ratio within the mass 
to charge ratio range of this calibration. This 
standard deviation may be combined with other 
statistical information to give an overall estimated 

10 error in a subsequent measurement. 

Although the method of calibration described above 
uses a least squares n th order polynomial, other 
calibration schemes may be employed such as, for 
example, non- linear and linear regression and splines. 

15 These may be more suited to particular types of data or 
mass analysers. In each case it is possible to extract 
information about the predicted error at particular mass 
to charge ratio values based on the statistics of the 
individual calibrant peaks. 

20 There are many potential sources of systematic 

error involved with mass measurement. However, many of 
these biases are constant during an experiment and can 
be removed from the final result by careful 
characterization of the system. This characterization, 

25 or calibration, is essential to ensure that the 

subsequent measurement is as accurate as possible. If a 
systematic error is identified that cannot be quantified 
then no reliable estimate of overall accuracy can be 
presented. 

3 0 Some sources of systematic error may not be 

constant during the course of a measurement or may not 
be accurately predicted and corrected for during 
calibration. An example of this is an unresolved mass 
interference . 
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Fig. 6 shows a mass spectrum wherein the 
statistical error is assumed to be negligible. Two mass 
peaks are shown corresponding to two different ion 
species. The two mass peaks are assumed to have a 
5 Gaussian distribution. The first mass peak PI has a 

mean mass to charge ratio value of Ml. The second mass 
peak P2 has a mean mass to charge ratio value of M2 . 
The intensity of the second peak P2 is set to be half 
that of the first mass peak PI. Both mass peaks P1,P2 

10 have a width at half height 5M of 0.05 mass units. The 
mass resolution is therefore 10,000 (FWHM definition). 
In the particular example shown in Fig. 6 Ml is 500.0000 
and M2 500.1500 i.e. the separation between Ml and M2 is 
0.15 mass units. It is clear that a resolution of 

15 10,000 is more than sufficient to separate the two 
different peaks P1,P2 and to enable accurate 
determination of the mean value of the peaks Ml , M2 . 

Fig. 7 illustrates a scenario wherein the 
separation between the two different mass peaks P1,P2 

20 has now been reduced from 0.15 mass units to only 0.05 
mass units e.g. Ml is 500.0000 but M2 is now 500.0500. 
It is apparent that the two mass peaks P1,P2 can not now 
be resolved from one another. The resultant single 
asymmetrical mass peak shown in bold has a mean mass to 

25 charge ratio value M r of 500.0167. The first mass peak 
PI may be an analyte ion which is intended to be 
measured whereas the second mass peak P2 may be an 
unexpected unresolved background or interference ion. 
Accordingly, the measurement of the mass to charge of 

3 0 the first (analyte) mass peak PI has an error E M i where: 



E = M r -Ml 
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In the particular example shown and described in 
relation to Fig. 7, the error E is 0.0167 da (33.4 ppm) . 

Mass interference effects causing systematic error 
can be recognized and in some cases approximately 
5 quantified by examination of the mass spectral peak 
shape compared to an idealized model of the expected 
peak shape at a particular mass to charge ratio value. 
This ideal peak shape may be characterized from that of 
standard reference peaks prior to a mass measurement. 

10 In some cases a Gaussian peak shape can be assumed. If 
the peak shape or width deviates beyond a pre -set 
criteria then the accuracy of the mass measurement of 
this mass to charge ratio value will be poor. The 
extent to which this effect can be accurately quantified 

15 will depend on the sophistication of the peak detection 
algorithm employed. 

Another source of systematic error can result when 
the signal intensity recorded exceeds the operational 
range of a signal amplification and/or a detection 

20 device. In certain mass analysers mass measurement 

accuracy can also be effected by space charge repulsion 
within the mass analyser itself. These effects can lead 
to shifts in the recorded mass to charge due to # for 
example, distortions in mass spectral peak shape as the 

25 intensity increases. The characteristics of the 

particular mass analyser and detection electronics 
utilized should ideally be characterized prior to 
measurement and the upper limit of the range at which an 
appropriate measurement of mass accuracy can be reported 

3 0 found . 

An additional source of systematic error can result 
from a change in the conditions under which the mass 
spectrometer is operating relative to those conditions 
existing when the system was calibrated. If a system 
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parameter has changed, resulting in the calibration 
becoming invalid, then the final mass measurement will 
be in error. This situation can occur, for instance, 
due to user intervention, the failure of a power supply, 
5 interference from an external source of electromagnetic 
radiation or due to mass drift caused by ambient 
temperature changes. The time period which has elapsed 
between calibration being performed or checked using a 
suitable reference compound and a mass measurement being 

10 reported are preferably monitored. The longer this time 
period becomes the less accuracy is preferably assigned 
to the final measurement . Although the error on the 
final measurement may not be readily quantifiable it is 
important that effort is made to recognize this 

15 situation and reflect a general uncertainty in accuracy 
with any measurement presented. 

To determine an overall figure of accuracy for an 
individual measurement, possible errors (described 
above) arising at each step in the mass measurement 

2 0 process may preferably be examined. 

Errors recorded in the instrument calibration 
procedure have been described above. Calibrations are 
generally made some time before a mass measurement of an 
analyte is performed. Due to the susceptibility of 
25 power supplies to temperature drift and the thermal 
expansion of components and materials used in the 
construction on mass analysers it is necessary to 
recognise and correct for this drift over time. 

To ensure that any minor drift in the instrument is 

3 0 accounted for when making a mass measurement an internal 

reference compound may be introduced at the same time or 
within a. short time period before or after a mass 
measurement is made. This not only allows the 
calibration expression to be adjusted for any small 
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changes but also provides an internal check of the 
stability of the system. For example, in a Time of 
Flight mass spectrometer a single point lock mass 
correction may be made. A reference material is 
5 introduced and the measured mass to charge ratio of a 
single mass spectral peak from this reference material 
is preferably monitored during the experiment. The 
empirical formula and hence the calculated mass to 
charge ratio for the internal calibrant is known. 

10 Changes in the overall gain of the calibration 

expression maybe made on a mass spectrum to mass 
spectrum basis. If more than one internal reference 
peak is present then more complex changes to the overall 
calibration expression may be incorporated including 

15 correction for any offset drift. 

An error in the determination of the mass to charge 
ratio value of the internal reference will result in a 
error in the mass measurement of the analyte . All the 
sources of error previously described should preferably 

2 0 be estimated for the measurement of the internal 
reference mass to charge ratio. 

In techniques where the secondary reference 
material is introduced at exactly the same time as the 
analyte compound, mass interference from background 

2 5 peaks or from analyte peaks may occur. In 

chromatographic techniques the analyte compound 
intensity will change with time as a chromatographic 
peak elutes. The internal reference compound, however, 
will remain at essentially constant intensity. Based on 

3 0 the intensity and the resolution of the internal 

reference peak the expected short-term scatter in the 
measured mass to charge ratio value for this may be 
calculated. If the mass to charge ratio value of the 
internal reference peak changes by significantly more 
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than the calculated short-term statistical scatter it is 
probable that a mass interference has occurred at the 
internal reference peak. This could be from mass 
interference due to the elution of the analyte compound 
5 or due to the elution of a contaminant compound 

unrelated to the analyte. Alternatively, the shift may 
be due to an instrument parameter change, failure of a 
power supply, or interference from an external source of 
electromagnetic radiation. If this shift is greater 

10 than the statistical scatter of the internal reference 

measurement, then the resultant error in the accuracy of 
the analyte peak may be estimated and indicated. Effort 
must be made to differentiate between long term drift 
corrected by an internal reference and short term 

15 instability causing systematic measurement errors. This 
mass shift may characteristically occur in a similar 
time-scale to the chromatographic peak elution width. 

Preferably, all the sources of statistical and 
possible systematic error previously described should be 

20 estimated for the analyte peak measurement. 

In addition, using chromatographic separation 
techniques several mass measurements of the analyte peak 
are preferably made during the chromatographic elution 
timescale. These are subsequently averaged and 

25 presented as a final reported mass to charge ratio. It 
is expected that the individual measurements across the 
chromatographic peak will have a statistical variation 
based upon their individual intensities. If a set of 
mass measurements from a peak from an eluting compound 

3 0 has an exceptionally large spread of values (a) compared 
to what is expected, then it may indicate the presence 
of interference from co-eluting or partially co-eluting 
compounds. In such situations no reliable estimation of 
accuracy may be reported. 
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If all the calculated errors above are determined 
in common units, e.g. mass units or parts per million 
(ppm) in mass to charge ratio, then the individual 
errors may be added in quadrature to yield an overall 
5 estimate of the errors associated with a single mass 
measurement : 



CJj =V C7 l 2 +0" 2 2 +CT 3 2 + Cr /i 2 

10 where <7i to a n represent the standard deviations 

calculated for each contribution to statistical error 
and a T is the total calculated standard deviation on the 
measurement of the analyte ion mass to charge ratio, 
value . 

15 The estimated accuracy of mass to charge ratio 

measurement is preferably recorded and optionally 
presented with each mass to charge ratio measurement 
displayed or listed in a mass spectrum. The figure 
reported preferably represents a scatter in multiples of 

20 ± <r T about the mean value. Based upon a Gaussian 

distribution the confidence in the measurement presented 
may expressed for different multiples of ± a T . The 
table below gives some common confidence limits derived 
from a Gaussian distribution. 

25 



Range 
Presented 


Probability of 
measurement 
falling within 
the range set 


± 0 . 5a T 


38 .29% 




68 .26% 


± 1 . 5a T 


86.64% 
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±2 o T 


95.43% ! 


± 2 . 5a T 


98 .76% 


± 3<7 T 


99.73% 


±3.5 a T 


99 . 95% 



A figure of ± 2 . 5 <r T representing a 98.76% 
confidence level may, for example, be presented. This 
figure may appear as an annotation within a mass 
5 spectrum or mass spectrum list and/or may be 

incorporated into the output from elemental composition 
calculation results . 

The accuracy of the measurement may also be 
indicated by automatically setting the number of decimal 

10 places for mass to charge ratio annotation that are 

presented to reflect the calculated confidence limit. 

It is important not to indicate too high a level of 
accuracy especially when restricting the minimum ppm 
window used for calculation of proposed possible 

15 elemental compositions. Situations in which 

insufficient data is available to calculate a meaningful 
estimation of the overall error in mass measurement 
should be determined. In the preferred embodiment an 
estimation of calculated error will only preferably be 

2 0 displayed if an internal reference compound was present 
during the measurement or correction was made within an 
acceptable time frame. 

The mass calibration used during the measurement 
should be valid. The criteria defining validity of the 

25 calibration should be determined for an individual mass 
spectrometer. These criteria include changes to 
critical instrument parameters which could affect 
calibration or if an unacceptable time has elapsed since 
the calibration was performed. 
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Fig. 8 and 9 illustrate the preferred method as 
applied to the exact mass measurement of a mixture of 
Polyethylene glycols ionised. using Electospray 
ionisation (ESI) . The resulting ions were analysed 
5 using an orthogonal acceleration Time of Flight mass 
spectrometer. Fig. 8 shows the residual errors after 
calibration. The x axis represents the calculated mass 
to charge value of the calibrant peaks based upon their 
known empirical formula. The y axis represents the 

10 difference between the calibrated mass to charge values 
and the known reference mass to charge ratio values 
after calibration had been applied. The y axis is 
displayed in units of parts per million error (ppm) . 
The vertical lines represent +/- 2a c where cr c is the 

15 statistical uncertainty in the peak center determination 
in ppm. The error bars represent the calculated scatter 
which could be expected from the reported mass to charge 
ratio from repetitive measurements of each calibration 
point within a 95% confidence limit. The dotted lines 

2 0 represent predicted 95% confidence bounds for subsequent 
measurements within the mass to charge ratio range of 
the calibration. These bounds were calculated based on 
the variance and covariance calculated using a third 
order weighted least squares linear regression. 

2 5 Fig. 9 shows a. mass spectrum according to the 

preferred embodiment of the mixture after calibration 
using an expression derived from the calibration shown 
in Fig. 8. The ion having a mass to charge ratio of 
745.41977 was used as single point internal secondary 

3 0 reference correction and hence no error band is 

presented for this particular ion. 

The predicted errors for each mass to charge ratio 
measurement were combined to provide an overall 95% 
confidence range for each measurement. The mass to 
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charge values displayed in the mass spectrum shown in 
Fig. 9 are annotated with the predicated 95% confidence 
ranges and are also truncated to three decimal places to 
reflect the predicted precision of the measurements. 
5 Although the present invention has been described 

with reference to preferred embodiments , it will be 
understood by those skilled in the art that various 
changes in form and detail may be made without departing 
from the scope of the invention as set forth in the 
10 accompanying claims. 



